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Abstract 



D . Provided the enhancement in the pp spectrum in radiative decay J/V' — > 

' 7PP observed by the BES collaboration is due to an existence of a pp 

molecular state, we calculate its binding energy and lifetime in the linear 

^ ■ a model. We consider a possibility that the enhancement is due to a pp 

H I resonance which is in either S-wave or P-wave structure and compare 

our results with the data. 



1 Introduction 

Recently, the BES collaboration has observed a near-threshold enhancement in the pp mass 
spectrum in the radiative decay J/ip -^ ^pp JL,. A similar report about the enhancement in 
B^ -^ D*^pp and B^ — > ppK^ decays has been published by the Belle collaboration ^j . 

There have been various interpretations for the observed enhancement. The enhancement 
can be understood if the final state interaction between p and p is properly considered, as some 
authors suggested ^. He et al. propose a possible mechanism that the final state of ^pp comes 
from an intermediate state of 7 + G where G is a '" or 0"^^ glueball ^. Meanwhile in analog to 
ao(980) and /o(980) which are supposed to be molecular states of KK, it is tempted to assume 
that pp constitute a bound state with quantum number '" or O"*""*". 



Under the assumption, one needs to evaluate the corresponding binding energy and Hfetime, 
then compare the theoretical results with the data. In this work, we employ the linear a model. 
Historically, there has been dispute about the linear a model where the cj-meson stands as a 
realistic scalar meson jSj whereas in alternative scenarios, it is suggested that the contribution of 
0" can be attributed to two-pion exchange _6 . In fact, the difference between the linear a model 
and non-linear a model is whether the 0^^ cj- meson is a substantial object, it corresponds to the 
linear or non-linear realization of the chiral lagrangian T. Because the low-energy QCD which 
is the underlying theory of hadron physics is fully non-perturbative, all the nonperturbative 
QCD parameters in the theory are so far not strictly derivable and have to be extracted by 
the data fitting. Therefore, determination of these parameters is somehow model-dependent 
and phenomenological. It is believed that at least for the leading order, all models would be 
applicable, even though they look somewhat different. As we employ the linear a model which 
is simpler in calculations, we take all the coefficients by fitting data. 

In our earlier work , we used the linear a model to calculate the properties of deuteron, and 
by fitting data we not only determine the value of ma but also fix the corresponding parameters 
of the linear a model. In this work we will use the same model with the parameters obtained 
by fitting the deuteron data to carry out all calculations for the pp bound state. 

The present BES data do not finally decide if the resonance is an S-wave or P-wave bound 
state, but only indicate that the position of the S-wave peak is below the threshold 2mp whereas 
the peak of the P-wave is a bit above the threshold. In our model, since the effective potential 
for the S-wave is attractive except a repulsive core near r ^ 0, the binding energy must be 
negative, so that the calculated mass of the S-wave bound state is below the 2mp threshold. 
Whereas for the P-wave due to the angular momentum barrier which is non-zero and positive, 
the binding energy becomes positive and the total mass is greater than 2mp. 

To evaluate the widths for both S-wave and P-wave, we investigate the dissociation mech- 
anism of the pp bound state. Since the central value of the mass of the S-wave bound state 
is smaller than 2mp, it dissolves into a pp pair via its width tail where the available energy is 
sufficient to produce a free pp pair. To evaluate the total width of the bound state, we need 
to achieve the imaginary part of the potential which is induced by the absorptive part of the 
loops in the pp elastic scattering amplitude (see the text for the concerned Feynman diagrams 
and some details). Thus according to the traditional method .0, , we derive the real part of the 
potential which mainly comes from the tree-level scattering amplitude where t-channel mesons 
are exchanged, including a, vr, p and uj. For the S-wave bound state not only t-channel exchange, 
but also the s-channel annihilation contribute. Namely in the s-channel, rj and rj' are the inter- 
mediate mesons and they contribute a real part and an imaginary part to the effective potential, 
the s-channel contributions are proportional to a delta function 6{r) in the non-relativistic ap- 
proximation. Thus the eigenenergy becomes Ejie — ^t &^d the time-factor is exp^—iEj^t — -^t) 
and the F corresponds to the total width and Ejie — i 2' i^ '^ solution of the Schrodinger equation 
with a complex potential. 

For the P-wave, the binding energy is positive and the angular momentum barrier pre- 
vents dissociation of the bound state. It is noted that since V'(O) = for the P-wave, the 
imaginary part of the complex potential which is proportional to 6{r), does not result in an 
imaginary part to the eigenenergy. The dissolution mechanism of the bound state is the quan- 
tum tunnelling. By the WKB approximation method ^j, the tunnelling transition proba- 



bility is exp[—2j^ ^/2^{V — E)dr], thus the total width of the P-wave bound state would be 
2fiexp[—2 J^ ^2\i{y — £^)dr], where // = ^ is the reduced mass. 

Substituting the potential no matter real and complex, into the Schrodinger equation, and 
solving it, one obtain both the eigenenergies and eigenfunctions of both S-wave and P-wave 
bound states. Then we can evaluate the masses and total widths of the bound states. That is 
the strategy of this work. 

This paper is organized as follows. After this introduction, we derive the formulation for 
the complex potential with a brief introduction of the linear a model. In sec. Ill, we substitute 
the potential into the Schrodinger equation and solve it to obtain the numerical result of the 
eigenenergy and eigenfunction, then obtain the masses and widths of the S- and P-wave bound 
states. In the section we also present all relevant parameters. The last section is devoted to our 
conclusion and discussion. 



2 The formulation 

(1) The necessary information about the model. 

In the linear a model, the effective Lagrangian is 

L = gip{a + -f5T ■7T)ip, (1) 

where ■0 is the wavefunction of the nucleon. When we calculate the scattering amplitude, we 
introduce a form factor to compensate the off-shell effects of the exchanged mesons. At each 
vertex, the form factor is written as [Hj 

' 7 , (2) 

A2 - g2 ' 

where A is a phenomenological parameter and its value is near 1 GeV. It is observed that as 
g^ ^ it becomes a constant and if A >> Mm, it turns to be unity. In the case, as the distance is 
infinitely large, the vertex looks like a perfect point, so the form factor is simply 1 or a constant. 
Whereas, as g^ — > oo, the form factor approaches to zero, namely, in this situation, the distance 
becomes very small, the inner structure (quark, gluon degrees of freedom) would manifest itself 
and the whole picture of hadron interaction is no longer valid, so the form factor is zero which 
cuts off the end effects. 

To derive an effective potential, one sets qq = and writes down the elastic scattering 
amplitude in the momentum space and then carries out a Fourier transformation turning the 
amplitude into an effective potential in the configuration space. Following the standard pro- 
cedure jH], we derive the effective potential from the scattering amplitude. Below, we present 
some details about the individual parts of the potential. 

(2) The effective potentials. 

(i) The real part of the potential. 



Here we first consider the meson exchanges at the t-channel, because the intermediate meson 
is a space-hke, it cannot be on its mass-shell, so that does not contribute to the imaginary part 
of the effective potential. Then we will go on discussing the s-channel contributions in subsection 
(ii). 

a. Via exchanging vr-meson: 

The effective vertex is 

L = gi^i^T ■ 7v^, (3) 

and obviously only vr" can be exchanged in our case. 
The scattering amplitude in the momentum space is 
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^-(q) = -T-^rT'^^32T(^i-q)(^2-q)(T^-f ) . (4) 



Following the standard procedure, we carry out a Fourier transformation on V^(q) and obtain 
the effective potential in the configuration space: 

K(r) = ^(^i-V)(^2-V)/.(r) (5) 

where 

-m^r -Ar / 2 _ \2\ -Ar 

Mr) = V-^+ ''8,A ^ («> 

b. Via a and p and to exchanges. 
The effective vertices are respectively 

La = gi^ixy, (7) 

Lp = g^^pht^r^'M'"', a = 1,2, 3, (8) 

Lu, = 5^^^rc.<V'7MV"^^• (9) 

The scattering amplitude via exchanging cj-meson is 

VM) = ^ 2f¥: 2, [4^^-V-q^-4(p-q)-»^-(qXp)](^^^°^' 
4m (q + m^J \A +q 

through a Fourier transformation, the potential is 
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Via exchanging vector- meson p{on\y p^ contributes), the effective potential is 
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For exchanging an u; vector meson, the expression is similar to that in the p case, but has an 
opposite sign to the p contribution due to the G-parity [^ H2] j thus one only needs to replace 
the corresponding parameter values, such as the mass and coupling constant for p by that for uj 
and add a minus sign in front of all the terms of Vp{r). For saving space, we dismiss the concrete 
expression for lo exchange. 

c. The real part of the potential 

A synthesis of all the individual contributions derived above stands as the real part of the 
effective potential, namely the traditional part of the effective potential as 

^e//(r) = V^{r)+V^ir)+Vp{r)+VUr) 

= Vo{r) + VLs{r) + Vpet{r) + VT + Vss. 



In the expression the leading part of the potential is 
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The relativistic correction which in our later numerical computations is treated as a perturbation 
to the leading part, is 
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The tensor potential is 
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d. The case of pp is different from the deuteron where the constituents are p and n, namely 
there is a pp annihilation at the s-channel, which would contribute a delta function to the 
effective potential. 

lipp is in S-wave with quantum number I'-'J^^ = 0^(0 ^), in the s-channel only a 0^(0 ^) 
meson can be exchanged. Here we only consider the lowest-lying pseudoscalar mesons of ^ i] 
and rj'. Their contribution can be written as 
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where m is the invariant mass \/{pi +^2)"^ and pi , p2 are the four- momenta of the constituents 
p and p respectively, and it is very close to 2mp. For the contribution of tj' , one only needs to 
replace the corresponding parameters values. It is noted that these contributions still belong 
to the real part of the effective potential. Below, we will derive the imaginary contributions 
induced by the absorptive part of loops at s-channel. 
(ii) The imaginary part of the complex potential. 

The corresponding Feynman diagrams are shown in Figs.l and 2. Fig.l is the self-energy of 
r] and ij' which are off-shell and Fig. 2 is a box diagram. Obviously, the elastic scattering of pp 
is a strong-interaction process, so that parity, isospin etc. quantum numbers must be conserved 
and as long as the pp bound state is of the '" structure, only ij and ij' can be exchanged in 
the s-channel (we neglect higher-resonances). 

The concerned couplings are ^2 ^1 E] 

Lpp^ = -l^ad^P-d^P, (10) 

Lvvp = 9vvpe^''^''d^V,{P-dxV„), (11) 

here P stands as pseudoscalar mesons, such as vr, r] and rj' etc. and V denotes vector mesons, 
such as oj and p etc. 

The imaginary part of the potential is obtained in the following way. First, we calculate the 
absorptive part of the loops by the Cutkosky cutting rule in the momentum space 14_ and carry 
out a Fourier transformation turning it into an imaginary part of the complex potential. 

(i) The contribution induced by the self-energy of 77 and -q' . 

We have obtained 
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Expanding the expression with respect to |p| and keeping terms up to p^, we have the final 
expression as 
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where the coefficient 0{m,m^,mri,A^) of p^ is obtained by a tedious but straightforward calcu- 
lation, here for saving space we ignore the details. 

Finally we obtain the imaginary part of the potential as 



It is also noted that for the P-wave resonance, the wavefunction at origin is zero, i.e. ip{0) = 0, 
at the leading order there is no s-channel contribution to the effective potential, and neither the 
imaginary part in the final effective potential, since in our approximation, all of them are pro- 
portional to (5^(r). 

In the practical computation, a popular approximation ^S] for the delta function 
is adopted. 

3 Numerical results 

By solving the Schrodinger equation, we obtain the zero-th order engenenergy and wave- 
function, where the L-S coupling and tensor terms are taken as perturbations and the imaginary 
part of the complex potential is treated in two ways. In terms of the traditional method of 
Quantum Mechanics, we can calculate the corrections. 

For the S-wave resonance we obtain the binding energy and the total width as 



m,(GeV) 


0.47 


0.48 


0.49 


0.50 


0.51 


0.52 


0.53 


0.54 


0.55 


A(GeV) 


0.59 


0.60 


0.61 


0.62 


0.63 


0.64 


0.65 


0.66 


0.67 


Es{MeY) 


-18.38 


-17.23 


-16.19 


-15.30 


-14.47 


-13.67 


-13.03 


-12.42 


-11.88 


r,(MeV) 


33.60 


30.08 


26.94 


24.15 


21.67 


19.49 


17.57 


15.90 


14.45 



Table 1: the theoretical results for the S-wave with perturbative method 

The number listed in table. 1 are obtained in terms of the perturbation method. Namely, we 
take the imaginary part of the complex potential as a perturbation as well as the L-S coupling 
and tensor terms, 

A^=<Vo|AF|Vo>, 

where AH = AH real + iAHimag and ipo is the wavefunction of zero-th order. 

As we sandwich the imaginary part of the complex potential between ijjq, the expectation 
value is the imaginary part of the complex eigenenergy as 



iAE, 
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and 



iVimir) = iAHi, 



Thus the total eigenenergy is 

E = Eq + AEj-eal + iAEimag 
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mag- 



Eq + AEreal ~ ^"^- 



Instead, one can solve the Schrodinger equation with a complex potential which would be 
divided into two coupled differential equations. The coupled equations cannot, in general, be 
solved analytically, but only numerically. We obtain the complex eigenenergy by solving the 
equation group and the results are listed in table 2. 



m,(GeV) 


0.47 


0.48 


0.49 


0.50 


0.51 


0.52 


0.53 


0.54 


0.55 


A(GeV) 


0.59 


0.60 


0.61 


0.62 


0.63 


0.64 


0.65 


0.66 


0.67 


Es{MeY) 


-14.10 


-13.64 


-13.17 


-12.71 


-12.25 


-11.81 


-11.38 


-10.99 


-10.62 


r,(MeV) 


21.90 


19.48 


17.29 


15.29 


13.48 


11.87 


10.45 


9.19 


8.11 



Table 2: the calculated results for the S-wave with direct program calculating 



Comparing the results in table 1 and 2, we find that there are some deviations of no more 
than 30%, and generally, the numbers obtained in the perturbation method are a bit greater than 
them by directly solving the coupled equations, but qualitatively, the two sets are consistent. 
Considering the experimental errors and theoretical uncertainties, we would conclude that the 
two sets of numbers agree with each other. 

For the P-wave, by solving the Schrodinger equation, we obtain the eigenenergy, and then 
by the WKB method, we can estimate the dissociation rate, which turns out to be the width of 
the P-wave resonance. The results are shown in Table 2. 



m^(GeV) 


0.47 


0.48 


0.49 


0.50 


0.51 


0.52 


0.53 


0.54 


0.55 


A(xlO-i GeV) 


8.25 


8.51 


8.76 


9.02 


9.28 


9.54 


9.81 


10.08 


10.35 


Ep{MeY) 


0.35 


0.43 


0.54 


0.63 


0.72 


0.80 


0.87 


0.94 


1.08 


rp(MeV) 


9.11 


11.42 


15.00 


17.96 


20.22 


22.64 


24.26 


25.91 


29.98 



Table 3: The calculated values for P-wave 



For the theoretical calculations, we have employed the following parameters as inputs: m 
0.938(GeV); m^ = 0.138(GeV); rup = 0.77(GeV); m^ = 0.783(GeV); m^ = 0.547(GeV); m^, 
0.958(GeV); m/g = 0.98(GeV); niao = 0.98(GeV) [E]; gNNw = QnNu = 13.5; qnnp = GnNu 

3.25 dZl; %^ = 0.4; ^al = o.f 
1.72(GeV-i)l7W/o = -9-01(GeV-i 
g^pp = -16(GeV-i)|I31|2ni. 



Irirjo 



4.11(GeV~i); 
-6.80(GeV-i; 



% 



rjrj cr 



T. 



TTTI ao 



2.65(GeV-i); 7^^/„ = 
= -7.80(GeV-i) [H; 



4 Conclusion and discussion 

In this work, in terms of the linear a model we investigate the spectrum and total width of 
the possible pp bound states. We consider two possibilities that the observed enhancement is 
due to a pp bound state in S- or P-waves respectively. 
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All the parameters employed in the calculations were obtained by fitting the data of deuteron. 
With the very precise measurement on the binding energy of deuteron and more or less accurate 
estimate of the s-d mixing and charge radius, there is only a narrow window in the parameter 
space |S|. Namely there is almost not much free room to adjust them, and neither is a large 
range for changing our theoretical calculations, as long as the model is employed. Therefore the 
newly observed resonance, if it is experimentally confirmed, can also provide an opportunity to 
further testify the linear a model. 

We derive the effective potential between p and p, and for the S-wave structure, we simply 
substitute it into the schrodinger equation to obtain the binding energy. We have also calculated 
the absorptive part of the concerned loops by the Cutkosky cutting rule and it becomes the 
imaginary part of the potential. Differently from the deuteron case which is a bound state of pn, 
for the pp case, there exist s-channel processes (see the figures in the text), which can contribute 
a real part (the tree level meson exchange) to the potential and an imaginary part through the 
loop diagrams in the channel and both of them are proportional to 5{t) at the concerned non- 
relativistic approximation. In this work, we ignore the dispersive part of the loops because it 
depends on the renormalization scheme and only makes a correction to the leading contribution 
of the real part of the potential, but keep the absorptive part which is the only source of the 
imaginary part of the potential. 

When we solve the Schrodinger equation with a complex potential, we have taken certain 
approximations to simplify the calculations. Then we obtain the mass and total width of the 
S-wave resonance. 

It is also noted that due to the G-parity structure of the A^A^ and A^iV systems where A^ 
refers to nucleons, the potentials contributed by a and p are of the same sign for the A^A^ and 
NN systems, but the potential induced by vr'' and a; should have opposite signs for the two 
systems |211l22j . In fact, for the deuteron case, which is in the pn structure, the contribution of 
vr*^ is repulsive, and so is that from p and u>. But for the pp system, vr*^ and u> induce attractive 
potentials while the contribution induced by a and p remain unchanged. It can qualitatively 
explain why the the binding energy for pp (about —18 MeV) is more negative than that for 
deuteron (about -2.22 MeV). 

For the P-wave case, where the wave function of origin is zero, i.e. ^(0) = 0, one does not need 
to calculate the s-channel contribution. The angular momentum barrier prevents dissociation 
of the bound state, but the quantum tunnelling leads to a final dissolution of the bound state 
and this tunnelling rate determines the total width or lifetime of the P-wave bound state of pp. 
Fig. 3 shows the effective potential for S-wave and P-wave respectively. The repulsive part at 
the region of small r is due to the vector-meson exchange. One also notes that for the S-wave, 
besides the repulsive core for small r the potential is attractive, whereas for the P-wave, there 
exists an angular barrier which results in an positive binding energy, i.e. the total mass of the 
P-wave resonance is above the threshold of 2mp. Since the barrier is not high, the binding energy 
is not far above zero and the total mass is very close to 2mp. To evaluate the tunnelling rate we 
use the WKB method, however, since the barrier is not much higher than the binding energy 
level, using the WKB method might bring up certain errors. Therefore the estimated width can 
only be valid to its order of magnitude. Indeed, with present experimental accuracy, we can 
satisfy ourselves with such numbers, but definitely the future experiments can provide us with 
much more information and by them we will modify our model and determine the concerned 
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parameters to higher accuracy. 

The newly observed enhancement by BES and Behe may have various interpretations, one of 
them is due to a resonance of pp. In this work we discuss this possibihty in the hnear o"— model 
and the obtained values are quantitatively consistent with the data. Our numerical results show, 
the total width and position of the proposed bound state, no matter S-wave or P-wave do not 
contradict the data, therefore both of them may be possible states which can accommodate the 
observed enhancement. 

The authors of P] suggested an alternative explanation, i.e. the final state interaction results 
in the observed enhancement. To decide which mechanism is right or dominant would wait for 
the future experiments. We hope that studies on the new resonance can enrich our knowledge 
about the hadron physics and the interactions at the hadron level. Our conclusion is that to 
confirm the observed enhancement, more precise measurements are needed. 
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